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Reheating in supersymmetric high scale inflation 
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Motivated by Refs we analyze how the inflaton decay reheats the Universe within supersym- 

metry. In a non-supersymmetric case the inflaton usually decays via preheating unless its couplings 
to other fields are very small. Naively one would expect that supersymmetry enhances bosonic pre- 
heating as it introduces new scalars such as squarks and sleptons. On the contrary, we point out that 
preheating is unlikely within supersymmetry. The reason is that flat directions in the scalar poten- 
tial, classified by gauge invariant combinations of slepton and squark fields, are generically displaced 
towards a large vacuum expectation value (VEV) in the early Universe. They induce supersymmetry 
preserving masses to the inflaton decay products through the Standard Model Yukawa couplings, 
which kinematically blocks preheating for VEVs > 10 13 GeV. The decay will become allowed only 
after the flat directions start oscillating, and once the flat direction VEV is sufficiently redshifted. 
For models with weak scale supersymmetry, this generically happens at a Hubble expansion rate: 
H ~ (10 -3 — 10 _1 ) TeV, at which time the inflaton decays in the perturbative regime. This is 
to our knowledge first analysis where the inflaton decay to the Standard Model particles is treated 
properly within supersymmetry. There are number of important consequences: no overproduction 
of dangerous supersymmetric relics (particularly gravitinos), no resonant excitation of superheavy 
dark matter, and no non-thermal leptogenesis through non-perturbative creation of the right-handed 
(s)neutrinos. Finally supersymmetric flat directions can even spoil hybrid inflation all together by 
not allowing the auxiliary field become tachyonic. 



I. INTRODUCTION 

Reheating after inflation connects the observable sector 
to the inflaton sector, which does not apriori carry Stan- 
dard Model (SM) charges. Often inflationary paradigm is 
realized with a SM gauge singlet inflaton whose origin and 
couplings to the matter cannot be explained within SM or 
its minimal extensions Therefore transferring infla- 
ton energy density into the SM degrees of freedom is the 
single most important phenomenon which would guaran- 
tee a successful big bang nucleosynthesis (BBN) [J]. 

The inflaton decay is the first and the most relevant 
part of reheating. Only one-particle decay of the non- 
relativistic inflaton quanta were considered initially Q. 
The treatment is valid if the energy transfer to the fields 
which are coupled to the inflaton takes place over many 
inflaton oscillations. This requires that the inflaton cou- 
plings to the SM fields are sufficiently small, which is 
justifiable if the inflaton (being a gauge singlet) couples 
through non-renormalizable operators to the matter sec- 
tor. Usually it is assumed that the plasma reaches com- 
plete kinetic and chemical equilibrium immediately after 
all the inflaton quanta have decayed, which is true in a 
non-supersymmetric scenarios @, H 0| ■ 

It was also pointed out that the coherent oscillations of 
the inflaton can create particles non-perturbatively 0, Q . 
This mechanism is called preheating and it is particularly 
efficient when the final products are bosonic degrees of 
freedom. It only takes about two dozens of oscillations to 
transfer the energy from the homogeneous condensate to 
non-zero modes of the final state(s) [9(. However, despite 
efficiently transferring the energy, preheating does not re- 
sult in a complete decay of the inflaton. In some cases the 



inflaton condensate fragments to form non-topological 
solitons which decays through surface evaporation [10l | . 
Irrespective of theses situations it is very challenging to 
understand thcrmalization of the preheated plasma. 

Although an epoch of perturbative reheating is an es- 
sential ingredient of any potentially realistic cosmolog- 
ical model [Til ] . preheating remains a possibility which 
can give rise to rich physical phenomena ranging from 
non-thermal production of particl es M , Q , production of 
gravitinos O^EEH] and moduli fl5||7 large isocurvature 
perturbations jl6| , amplifyin g gr avity waves [l7j ] , genera- 
tion of large non-Gaussianity [18L [1911 . non-thermal source 
for leptogenesis and baryogenesis [20( , and the formation 
of topological and non-topological solitons [13, HH • 

Supersymmetry (SUSY) provides the framework for 
most widely studied extensions of the SM physics be- 
yond the electroweak scale. Therefore, given the large 
energy density of the inflaton, it is then pertinent to ask 
how the inflaton decays within SUSY. Particularly, into 
the relativistic degrees of freedom of the minimal super- 
symmetric SM (MSSM). Then the next question is how 
quickly do they thermalizc. This very last question was 
recently addressed in Ref. Q. 

In a Non-SUSY case it is known that the inflaton decay 
products thermalize very quickly because of the efficiency 
of interactions mediated by the massless gauge bosons of 
the SM ;6]. Therefore the reheat temperature is mainly 
governed by the inflaton decay width. In Refs. p| [If, 
we pointed out that the process of thermalization is in 
general painstakingly slow within SUSY. The Universe 
loiters in a phase of a quasi-thermal equilibrium after 
the decay of the inflaton, which results in a very low 
reheat temperature, i.e. T R ~ O(10 3 - 10 7 ) GeV. The 
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final reheat temperature does not depend on the decay 
width of the inflaton rather on a thermalization time scale 
which depends on the vacuum expectation value (VEV) 
of the squarks and sleptons. 

What was not addressed in Ref. Q is the initial phase 
of the inflaton decay. In this paper we will fill up that 
gap. We would like to know how the inflaton decays: 
whether perturbatively or non perturbatively. This simple 
question is so relevant that depending on the nature of 
the inflaton decay there would be different consequences 
all together. For instance production of baryons, cold 
dark matter, magneto-genesis, electroweak baryogenesis, 
production of dangerous relics and their abundances de- 
pend on the nature of the primordial plasma. Therefore 
this is an important topic which relates the early Uni- 
verse physics to phenomenology. 

An important fact is the presence of flat directions 
along which the scalar potential identically vanishes in 
the limit of exact SUSY. In the MSSM alone there are 
nearly 300 flat directions [13], which are made up of 
gauge invariant combinations of squarks, sleptons and 
Higgses. These are none but the simplest examples of 
moduli near points of enhanced symmetry. During infla- 
tion more than one MSSM flat directions (orthogonal in 
flavor basis) 12311 are expected to develop large VEVS, for 
a review, see |24| . 

A large VEV of the MSSM flat directions, during 
and after inflation, spontaneously breaks the SM gauge 
group, and gives masses to the gauge bosons and gaug- 
inos similar to the Higgs mechanism. Many of the flat 
directions break the entire SM gauge group 25]. The 
flat direction VEV also induces large SUSY preserving 
masses to (s)quarks, (s)leptons and Higgs/Higgsino fields 
during and after inflation. As we will see, such large 
masses kinematically prohibit non-perturbative inflaton 
decay into MSSM fields. The initial stage of inflaton os- 
cillations thus produces no significant fraction of MSSM 
particles Q. The decay will occur much later after the 
fiat direction oscillations have started and their VEV has 
been sufficiently redshifted. Such a delayed decay of the 
inflaton is typically in the perturbative regime. 

Rest of the paper is organized as follows. In Section 
II we briefly discuss inflaton decay in a non-SUSY case. 
We highlight inflaton couplings to the MSSM fields in 
Section III, and MSSM flat direction couplings to the 
inflaton decay products in Section IV. We then discuss 
the flat direction dynamics during and after inflation in 
Section V. In Section VI we explain how the VEV of flat 
directions prevent non-perturbative decay of the inflaton. 
We illustrate in Section VII how the inflaton eventually 
decays perturbatively, and discuss various cosmological 
consequences in Section VIII. Finally we briefly mention 
our conclusion. We have added appendices discussing 
some minute details for the paper to be self-contained. 



II. INFLATON DECAY IN A NON-SUSY CASE 

First let us briefly review the initial stage of the infla- 
ton decay which is typically non-perturbative, i.e. pre- 
heating, in a non-SUSY case. Our focus is on bosonic 
preheating which acts most efficiently in transferring the 
energy density from the inflaton oscillations. We consider 
models of large field inflation, such as chaotic inflation, 
for which bosonic preheating is most pronounced. The 
relevant renormalizable couplings between the inflaton <f> 
and a scalar field x will read from the following potential: 

V = + a^x 2 + h 2 0\ 2 + A X 4 , (1) 

where we have considered cj) and x to be real. Here a 
is a coupling which has a [mass] dimension. The only 
scalar field in the SM is the Higgs doublet. Therefore in 
a realistic case x denotes the real and imaginary parts 
of the Higgs components 1 . Note that the cubic interac- 
tion term is required for a complete inflaton decay. The 
quartic self-coupling of x is required to bound the poten- 
tial from below along the x direction. The dimensionless 
couplings <j/m,p and h (as well as A) are not related to 
each other, hence either of the cubic or the quartic terms 
can dominate at the beginning of inflaton oscillations (i.e. 
when the Hubble expansion rate is H(t) ~ and the 
amplitude of oscillations is <f> ~ O(Mp)). 

• a <C h 2 Mp: In this regime the h 2 cf> 2 x 2 term is dom- 
inant at the beginning of the inflaton oscillations. 
This case has been studied in detail in first two ref- 
erences of @ . For a nominal value of the inflaton 
mass, = 10 13 GeV, non-perturbative x produc- 

/ -V /2 

tion with a physical momentum, k <J I hm^cf) 1 , 

takes place if h > 10~ 6 . Particle production is par- 
ticularly efficient if h > 3 x 10~ 4 , and results in 
an explosive transfer of energy to x quanta which 
ends when re-scatterings destroy the inflaton con- 
densate 2 . The whole process happens over a time 
scale ~ lSOm^ 1 , which depends logarithmically on 
h 3 . 



1 Since the SM fermions are chiral, the inflaton can only couple to 
them through dimension-5 operators. The same holds for cou- 
pling to gauge bosons where the inflaton is coupled to gauge 
field strengths. Such couplings are non-renormalizable and sup- 
pressed compared to those in Eq. JT}, and hence negligible (we 
consider the scale of non-renormalizable operator is governed by 
M P ). 

2 Further note that the h 2 (f> 2 x 2 term does not produce any signif- 
icant non-Gaussianity if h > 10 — 5 , see the first reference of [ljjj . 

3 In a non-SUSY case efficient preheating happens over a narrow 
window 3 X 10~ 4 < h < 10~ 3 . The reason is that the h 2 (j> 2 x 2 
term yields a quartic self-coupling for the inflaton at a one-loop 
level which is constrained by the COBE normalization of the 
density perturbations [53. l23j. However in SUSY this correction 
is canceled out by that from fermionic partner of x> so in princi- 
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• a ^> ft 2 Alp: In this regime the cubic term <j4>x 2 
dominates. This case was recently considered in 
Refs. [3(3, [3l|, where the the x field becomes tachy- 
onic during half of each oscillation. For a > 
m^/Mp (which amounts to a > 10 7 GeV for 

rritf, = 10 13 GeV) this tachyonic instability trans- 
fers energy from the oscillating condensate very ef- 
ficiently to the \ quanta with a physical momentum 

( V /2 

k <J I a<f> \ . Particle production ceases when the 
back-reaction from x self-coupling induces a mass- 
squared ;> crip. Depending on the size of A, most of 
the energy density may or may not be in x quanta 
by the time backreaction becomes important [31[ . 

Couple of points to note here. In the borderline regime 
a ~ h 2 Mp, the cubic and quartic interaction terms 
are comparable. The inflaton decay happens due to a 
combination of resonant and tachyonic instabilities. If 
h -C m^/Mp and a <C m^/Mp, the inflaton decays per- 
turbatively via the cubic interaction term. However this 
requires very small couplings: h, (cr/m^) < 10~ 6 . There- 
fore, unless the inflaton is only gravitationally coupled to 
other fields, the initial stage of its decay will be generi- 
cally non-perturbative. 

The plasma from the non-perturbative inflaton decay 
eventually reaches full thermal equilibrium, though, at 
time scales much longer than that of preheating itself [26|, 
l27l . [28[. The occupation number of particles is fk ^> 
1 in the meantime. This implies that dangerous relics 
(such as gravitino and moduli) can be produced much 
more copiously in the aftermath of preheating than in full 
thermal equilibrium P, [l5|, [28| . This is a negative aspect 
of an initial stage of preheating. One usually seeks a 
late stage of entropy release, in order to dilute the excess 
of relics. As we shall show, supersymmetry naturally 
provides us a tool to undo preheating completely. 



III. INFLATON COUPLINGS TO MATTER IN 
SUSY 

Inflaton couplings to (MS)SM fields is of utmost im- 
portance for (p)reheating. In all the relevant papers, for 
instance see Refs. [1, Q, inflaton couplings to matter has 
not been dealt with carefully. Only toy models have been 
considered which have no relevance to SM physics 5 . The 



pie one could expect a rather broader range of parameter space 
within SUSY. 

4 In Ref. [3ll . the authors attempted to motivate the cubic cou- 
pling from SUSY. However they missed vital ingredients which 
exist in a realistic case, such as the strength of the SM couplings, 
contributions from SUSY flat directions, etc. In this paper we 
wish to note that neither the couplings nor the interactions are 
taken arbitrarily. 

5 Fermionic preheating has been discussed in Refs. |44fl , but the 
importance of SM gauge invariance was grossly neglected. In fact 



importance of gauge invariance was first highlighted in 
Refs. 0,0, SIP 6 . 

In almost all known F and/or D-term models of infla- 
tion the inflaton, 0, is considered to be an absolute gauge 
singlet. Then the main question arises how the inflaton 
couples to the matter. This is one of the most pertinent 
issues which connects inflation to a hot big bang cosmol- 
ogy 

First note the field content of MSSM which is governed 
by the following superpotential: 

Wmssm = A u QH u u + A d QH d d + A e LH d e + fxH. u H d , 

(2) 

where H u , H^, Q, L, u, d, e in Eq. ([2]) are chiral super- 
fields representing the two Higgs fields (and their Hig- 
gsino partners), left-handed (LH) (s)quark doublets, 
right-handed (RH) up- and down-type (s)quarks, LH 
(s)lepton doublets and RH (s)leptons respectively. The 
dimensionless Yukawa couplings A u , A^, A e are 3x3 ma- 
trices in the flavor space, and we have omitted the gauge 
and flavor indices. The last term is the [i term, which is 
a supersymmetric version of the SM Higgs boson mass. 

There exist two gauge-invariant combinations of only 
two superfields: 

H u Hd , H U L. (3) 

The combinations which include three superfields are: 

H u Qu , H d Qd , H d Le , QLd , udd , LLe. (4) 

SUSY together with gauge symmetry requires that the 
inflaton superfield be coupled to these combinations 7 . 
The terms 3>H u H c i and <I>H U L have dimension four, and 
hence are renormalizable. On the other hand, the interac- 
tion terms that couple the inflaton to the combinations in 
Eq. ((4]) have dimension five and are non-renormalizablc. 
In following we focus on renormalizable interactions of 
the inflaton with matter which play the dominant role in 
its decay 8 . 



the inflaton couples to SM fermions through non-renormalizable 
dimension 5 operators, and therefore preheating into SM fermion 
is unlikely. 

6 In Refs. [34l | the inflaton belonged to a gauge sector which can 
carry SM charges based on the ideas of assisted inflation [35| . In 
which case the inflaton couplings to matter are governed by the 
usual Yukawas. However it is hard to construct realistic models 
without gauge singlets. 

7 It is possible that the inflaton mainly decays to another singlet 
(for example, the RH neutrinos) superfield, see the discussion in 
section IXI Bl 

8 We note that terms representing gauge-invariant coupling of the 
inflaton to the gauge fields and gauginos are also of dimension 
five, and hence preheating into them will be suppressed. 
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A. Two choices of renormalizable couplings 

The simplest case is when the inflaton is coupled to 
matter via superpotential terms of the form: 

2 5 *H u H d , 2 3 $H U L. (5) 

where g can be as large as 0(1). The factor of 2 
as we shall see, leads to convenience in field redefini- 
tions. Besides the SM gauge group the MSSM La- 
grangian is also invariant under a discrete Z2 symme- 
try namely "i?-parity" . This symmetry assigns the num- 
ber R = (— 1) " to the component fields where 

B, L, S denote the baryon number, lepton number and 
spin respectively. This amounts to +1 for the SM fields 
and — 1 for their supersymmetric partners. As a result 
the lightest supersymmetric particle (LSP) will be stable 
and can account for dark matter in the Universe. This 
is one of the most remarkable cosmological features of 
MSSM. 



B. Preserving imparity 



irrelevant during inflaton oscillations. We have neglected 
the inflaton coupling to the fermionic partners of x as 
we focus on the bosonic preheating here. However our 
analysis will follow similarly to the fermionic case and 
the same conclusions hold for fermionic preheating too. 
In addition to the terms in Eq. ([9]) there are also 

the self- and-cross-couplings, (g 2 /4) (xl ~ xl) + &xlxl, 
arising from the superpotential and D-terms respectively 
(a is a gauge fine structure constant). Therefore even in 
the simplest SUSY set up the scalar potential is more 
involved than the non-SUSY case given in Eq. (JXJ) , which 
can alter the picture of preheating presented in the lit- 
erature [37], ■ Note however that these terms become 
important after particle production has started. Here we 
focus on the terms in Eq. ([9]) which are relevant for par- 
ticle creation from the very beginning of the oscillations. 

A remarkable feature in Eq. ([9]) is that SUSY naturally 
relates the strength of cubic 4>x 2 and quartic (j> 2 x 2 inter- 
actions. We re-emphasize that the cubic term is required 
for complete decay of the inflaton field. This is a natural 
consequence of SUSY which holds so long as the inflaton 
mass is larger than the soft SUSY breaking masses. 



Preserving i?-parity at the renormalizable level further 
constrains inflaton couplings to matter. Note that H u H f ; 
is assigned +1 under i?-parity, while H„L has the oppo- 
site assignment — 1. Therefore only one of the couplings 
in Eq. ([5]) preserves i?-parity: ^fL u lld if i?# — +1, and 
3>H„L if i?# = — 1 (such as models where the RH sneu- 
trino plays the role of the inflaton (36|). Therefore the 
renormalizable inflaton coupling to matter can be repre- 
sented as 



where 



2g*H u * 



* = H n if i? # = +1 , 
* = L if = -1 



(6) 



(7) 



Taking into account of the inflaton superpotential mass 
term: and after defining 



Xi. 



(Hu ± g) 

V2 



(8) 



and with the help of Eqs. (|24I25[) . see appendix IXI A[ 
we find the renormalizable part of the potential which is 
relevant for the inflaton decay into MSSM scalars is given 
by: 



g 2 (t> 2 x 2 ± ~/|fl' m 



(9) 



where \ denotes the scalar component of X12 superfields, 
and we have only considered the real parts of the inflaton, 
0, and x field. Further note that the cubic interaction 
term appears with different signs for xi and X2, but this is 



IV. FLAT DIRECTION COUPLINGS TO 
INFLATON DECAY PRODUCTS 

Consider a MSSM flat direction, ip, with the corre- 
sponding superfield denoted by <p (only for flat directions 
we are denoting the superfield and the field with the same 
notation). For a brief discussion on MSSM flat directions, 
see Appendices IXI AllXI CI Note that the ip and X super- 
fields are linear combinations of the MSSM superfields, 
see Eq. ([8]), and hence are coupled through the MSSM 
superpotential in Eq. (|2|) . The couplings are nothing but 
the (MS)SM Yukawas. Then the MSSM superpotential 
can be recast in the following form: 



(10) 



where £1.2 are some MSSM superfields 9 . For example 
consider the case where ip is a flat direction classified by 
the udd monomial, and ^ = H c ;. In this case Si, 2 are Q 
superfields and Ai,2 correspond to X u and Xd respectively, 
see Eq. ©. After using Eq. © we find: 



W D ^X^Ei + ^LX^ 
y/2 \/2 



This results in: 



uda 2 mV , A 



2\ 1/2 



(11) 



(12) 



where we have again considered the real part of x- 



Note that Si 7^ ^ and S2 7^ H u , since ip is a non-gauge-singlet. 
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Let us determine the strength of flat direction cou- 
pling to X: which is denoted in Eq. (fT2f by A. Note that 
the first generation of (s)leptons and (s)quarks have a 
Yukawa coupling ~ O(10~ 6 — 10~ 5 ), while the rest of 
the SM Yukawa couplings are > 1(T 3 . The MSSM flat 
directions can be grouped in 6 categories mentioned in 
Appendix IXI C| out of which 

• Only 11 directions: 3 udds, 6 QdLs, 1 LLddd and 

1 LLe have couplings to MSSM particles/sparticles 
such that A < 0(1CT 5 ), 

• Rest of the flat directions have A > 3 x 1CP 4 . 



V. FLAT DIRECTION POTENTIAL 

The flat directions are lifted by soft SUSY break- 
ing mass term, m ~ O(TeV), Hubble induced cor- 
rections and superpotential corrections of type: W ~ 
X n ip n /nM n ~ 3 [H 

I |2(n-l) 

V D (m 2 + c H H 2 ) \p\ 2 + \V{^—^ , (13) 

with n > 4. Here M is the scale of new physics which in- 
duces the non-renormalizable terms, typically the Planck 
scale M = Mp or the grand unification scale M — Mqut ■ 

Note that ch can have either sign. If ch J> 1, the flat 
direction mass is > if. It therefore settles at the origin 
during inflation and remains there 10 . Since (ip) = at 
all times, the flat direction will have no interesting con- 
sequences in this case. However there is a large class of 
theories which predicts ch < and also cr <C +1. Neg- 
ative ch may arise naturally if the inflaton and MSSM 
flat directions have positive higher order couplings in 
the Kahler potential, i.e. $t$^y; [39], such that all 
the eigenvalues of the Kahler matrix is positive definite. 
There is no symmetry which prohibits such couplings. 

Moreover string theory, which we believe will provide 
the true low energy effective theory, also generically pre- 
dicts no-scale type Kahler potential based on Heisenberg 
symmetry [40| . which at tree level gives no correction to 
the flat direction mass, i.e. c# = 0. However a Hub- 
ble induced mass term is generated at a loop level, be- 
cause MSSM superpotential (as well as D-terms) breaks 
the Heisenberg symmetry, which induces a calculable but 
small contribution, i.e. ch < 10 -2 [4l[. Moreover, even 
starting with ch > at a high scale, it is possible that ch 
quickly changes sign due to loop corrections from large 
Hubble-induced SUSY breaking terms [43[ . 

In the absence of c# J> 1, the flat direction remains 
flat during inflation as the Hubble expansion rate is 
Hi 3> ttiq . Therefore quantum fluctuations are free to 



This has a similar origin as a supergravity inflationary r;-problcm, 
see 0. 



accumulate (in a coherent state) along ip and form a 
condensate with a large VEV, ipo- Because inflation 
smoothes out all gradients, only the homogeneous con- 
densate mode survives. However, the zero point fluctua- 
tions of the condensate impart a small, and in inflation- 
ary models a calculable, spectrum of perturbations on 
the condensate [24| . 

If the higher-order superpotential term is forbidden, 
due to an R— symmetry (or a set of R— symmetries) [42j . 
then we naturally have, <po ~ Mp [3{|. On the other 
hand, pa -C Mp will be possible if non-renormalizable 
superpotential terms are allowed. As shown in [221 ] . in the 
absence of any i?-symmetry, all the MSSM flat directions 
are lifted by higher-order terms with n < 9. If a flat 
direction is lifted at the superpotential level n, the VEV 
that it acquires during inflation cannot exceed: 

^-(ffxM"- 3 ) 17 ^ , (14) 

where Hi is the expansion rate of the Universe in the 
inflationary epoch. 

After inflation, H(t) cx i -1 , the flat direction stays 
at a relatively larger VEV due to large Hubble friction 
term, note that the Hubble expansion rate gradually 
decreases but it is still large compared to m . When 
H(t) ~ mo, the condensate along the flat direction starts 
oscillating around the origin with an initial amplitude 

ip- m ~ (mrjMp -3 ) 1 " 2 \ From then on (cp) is redshifted 
by the Hubble expansion cx H for matter dominated and 
cx ii 3 / 4 for radiation dominated Universe. 



VI. NO PREHEATING IN SUSY 

In order to understand the preheating dynamics it is 
important to take into account of x coupling to the infla- 
ton (j), as well as to the MSSM flat direction, tp, which is 
displaced away from its minimum (towards large VEVs) 
during inflation. The governing potential can be obtained 
from Eqs. (j9TT2"|) 




As mentioned in the previous section, we generically have 
A > 3 x 10~ 4 , and g can be as large as ~ O(l). 

After mode decomposition of the field \i the energy 
of the mode with momentum fc, denoted by Xki is given 
by 1: 

1 /2 

to k = (k 2 + 2g 2 {4>) 2 + V2g m(i> ((t>) + 2A 2 (<^) 2 ) . (16) 

We have frozen the expansion of the Universe. Including 
the expansion will not change our conclusions anyway. 
First note that during inflation the inflaton VEV is large, 
i.e. ((f) > Alp. Therefore if g > 1CP 6 the inflaton induces 
a large mass g(4>) > Hi for x during inflation. 

As a result, Xi quickly settles down to the minimum, 
i.e. (x) = 0, even if it is initially displaced, and remains 
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there. Therefore, <p, does not receive any mass correc- 
tions from its coupling to x during inflation. Note that 
the VEV of the flat direction, ip, induces a large mass, 
Ai^o, to the x field during inflation. 

At the end of inflation, i.e. when H{t) ~ m^, the in- 
flaton starts oscillating with frequency m$ and an initial 
amplitude 0(M P ). Note that for g > 10~ 6 the quartic 
inflaton coupling takes over the cubic one in Eq. (fTS)) . 
In the interval mo < H(t) < the flat direction 
VEV slides very slowly because of the under damped 
motion due to large Hubble friction term, the flat di- 
rection effectively slow rolls. Non-perturbative produc- 
tion of x quanta will occur if there is a non-adiabatic 
time- variation in the energy, i.e. that duik/dt J> w 2 . The 
inflaton oscillations result in a time- varying contribution 
to LUk , while the flat direction coupling to x yields a vir- 
tually constant piece. 

Obviously the piece induced by the flat direction VEV 
weakens the non-adiabaticity condition 11 . Indeed time- 
variation of ujk will be adiabatic at all times 



provided that 

A 2 (<^) 2 > #m , (18) 

where <j> is the amplitude of the inflaton oscillations. We 
find the most conservative bound by considering the most 
optimistic situation for preheating: 

• The largest possible amplitude for the inflaton os- 
cillations, (fi r*> 0{Mp). 

• The largest possible coupling to the inflaton, g ~ 
0(1). 

We remind that in the absence of flat direction VEV 
preheating would enter an explosive stage when H(t) ~ 
1(T 2 TO0, at which time $ ~ 1CT 2 M P @. Therefore, 
to prevent efficient preheating, it will be actually suf- 
ficient to satisfy Eq. (jTTJ) at this time rather than the 
very beginning of inflaton oscillations. Moreover, as men- 
tioned earlier, having g ~ 0(1) also implies a large x 
self-coupling in supersymmetry. Then one in addition 
expects preheating to be considerably suppressed due to 
self- interactions (33,[1^|. Nevertheless we want to find the 
strongest bound on the flat direction VEV which shuts-off 
non-perturbative particle production at all times and for 
the largest coupling to the inflaton. The natural conclu- 
sion is that there will be no preheating in more realistic 
situations. 



It turns out from Eqs. (|16ll7p that the energy of 
mode Xk changes adiabatically at all times if Eq. (fl"8]) 
is satisfied at the beginning of inflaton oscillations, i.e. 
H = rri(p 12 . Hence there will be no resonant production 
of x quanta, provided that 

<Po > A" 1 (<?M P m ) 1/2 . (19) 

This surmounts to a kinematical blocking of preheating 
by inducing a piece (which is virtually constant at time 
scales of interest) to the mass of inflaton decay products 
due to their couplings to a flat direction which has a large 
VEV. 

Once g<fi ~ m^, the Hubble expansion rate becomes: 
H(t) ~ rrA/gMp, and the cubic interaction term (fix 2 
takes over the quartic one. For = 10 13 GeV and 
g > 10~ 6 this happens when the expansion rate is still 
> 1 TeV. Eq. implies that \ 2 ipl > m|, and hence 
the flat direction VEV totally dominates the mass of x 
when the cubic term has taken over. Therefore there will 
be no tachyonic instability in the mass of x which would 
otherwise occur during half of each inflaton oscillation 
and could lead to efficient particle production [3l|. 

Now let us find the range of VEVs for which a MSSM 
flat direction can satisfy the condition given in Eq. (fT9"l) 
and shut-off preheating. We choose a nominal value of 
the inflaton mass — 10 13 GeV. As discussed in the 
previous subsection, we have A > 3 x 1(T 4 for all flat 
directions but few exceptional ones. The required values 
of (fo are depicted in Table. 1 for A 2 = 1CP 7 and A = 
l) -1 . The three categoric values of inflaton coupling g — 
1, 3 x 1CP 4 , 10~ 6 illustrate the distinctive regimes where 
inflaton decay would take place in the absence of flat 
direction VEV: 

• 3x10 <(?<1: In this case preheating via 
the quartic interaction term g 2 <fi 2 X 2 would be effi- 
cient @. 





5 = 1 


g = 3 x 1CT 4 


9 = 10- 6 


A 2 ~ lCr 7 


tpo ~ Mp 


ip > 3 x 10 17 


<A) > 10 16 


A 2 ~ lCr 1 


<po > 10 16 


tpo > 3 x 10 14 


<Po > 10 13 



Table 1: The flat direction VEV, inflaton coupling g, 
and two values of generic flat direction (Yukawa) cou- 
plings A within MSSM are illustrated. The VEVs are 
denoted in GeV. For the VEVs above the quoted num- 
bers preheating is kinematically blocked. 



11 Note that in an absence of 2\ 2 (ip) 2 term in Eq. ( 1161 . the adi- 
abaticity condition would be violated every time the inflaton 
would cross the origin leading to a copious production of x par- 
ticles with momentum k < (g(j>m$\ @, [t| (<j> is the amplitude 
of the inflaton oscillations) . 



In the absence of resonant particle production the amplitude 
of the inflaton oscillations is redshifted like <j> oc H(t) due to 
the Hubble expansion. While, for H > mo, we have (tp) oc 
fl-(t)l/(n-3) j see Eq JJU Hence the RH side of Eq. JTS) be- 
comes increasingly larger than its LH side, and the adiabaticity 
condition will be satisfied more comfortably as time goes by. 
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• 10 -6 < g < 3xlO~ 4 : In this case preheating via the 
quartic term would not be efficient, however, the 
cubic interaction term (jm^ / \/2) <ftx 2 could lead 
to efficient particle production [3l| . 

• g < 1CV 6 : In this case inflaton would decay in the 
perturbative regime from the beginning, thus no 
resonant and/or tachyonic particle production. 

It is important to note that preheating is always shut- 
off for a sub-Planckian fiat direction VEV. Even in the 
most extreme case, i.e. the largest inflaton coupling g ~ 
0(1) and for a flat direction coupling A ~ 1CV 4 , the 
required VEV for the flat direction is ipo — Mp. Much 
smaller VEVs, ipo > 10 13 GeV, are required for moderate 
values of g and/or A ~ 0(1) (which is the case for flat 
directions including a sizeable component of Q and/or u 
from the third generation). 

A. Additional observations 

Some words on fermionic preheating, see Refs. [44| . 
Preheating of superheavy fermions can be much more 
efficient than bosonic preheating (20| . However within 
SUSY this is not the case as the symmetry between 
bosons and fermions implies similar equations for the mo- 
mentum excitations, see Eq. (fT5|) . The same also holds 
for kinematical blocking of preheating by the flat direc- 
tion VEV as the induced mass \(<p), is SUSY preserving. 
Therefore, the condition for the adiabatic change in vac- 
uum remains as it is in the case of bosonic excitations, 
and fermionic preheating is prohibited for the flat direc- 
tion VEV given in Eq. gg) . 

Finally let us compare the flat direction VEV obtained 
in the presence of non-renormalizable terms with those 
given in Table. 1. For = 10 13 GeV, and M — M P in 
Eq. (fT5|) . we obtain from Eq. JTJJ the field values shown 
in table. 2 right at the end of inflation (i.e. at H ~ m$). 





n = 4 


n = 6 


n = 9 


Po 


10 16 GeV 


10 17 GeV 


10 18 GeV 



Table 2: The flat direction VEV at the beginning of 
inflaton oscillations in presence of the non-renormalizable 
term in the potential, see Eq. (fT3"j) . 



These values comfortably lie within the range de- 
picted in Table. 1. For H < the amplitude of in- 
flaton oscillations is redshifted <f> cx H, while the flat 
direction VEV slides down to an instantaneous value 
{if) ~ (H(t)M n - 3 ) 1/n ~ 2 , see Eq. (JHJ) . Therefore if 
A 2 c/?q > gMprritf, at the onset of inflaton oscillations, we 
will have A 2 (ip) 2 > g4>m^, at later times (note that n > 4). 
This implies that even when the flat direction VEV slowly 



slides down in a non-renormalizable potential preheating 
will remain shut-off. 

Once H(t) ~ uiq ~ 1 TeV, soft SUSY breaking mass 
term in the potential takes over and the flat direction 
starts oscillating around its origin with an initial ampli- 
tude: tp- m ~ (moA/p -3 ) 1 ^" 2 . However this happens at 
time scales hierarchically longer than those relevant for 
the preheating phenomena. 

VII. INFLATON LIFE TIME 

Since preheating is kinematically shut-off, the inflaton 
simply oscillates with a decreasing amplitude due to the 
Hubble expansion rate for H(t) < m$. The Universe is 
therefore dominated by the inflaton oscillations, which 
for m 2 ^4> 2 potential act as a non-relativistic matter, im- 
plying that <j) oc H(t). 

Once H(t) ~ mo, the flat direction also starts oscillat- 
ing and, due to Hubble damping, (ip) cx H. We remind 
that for g(f> <C m^, the cubic interaction term 4>x 2 is dom- 
inant and production of x can only occur in a narrow mo- 
mentum band peaked around k = [1, Q. Note that 
at the onset of flat direction oscillations <p = "'noMp/m^, 
and hence <?</> <SC m<j, even if g = 1. Therefore the inflaton 
decay will be kinematically forbidden until X(cp) < 

In the universe which is filled by the inflaton oscilla- 
tions (ip) cx H(t) for H < mo- Thus the inflaton eventu- 
ally decays when the Hubble expansion rate becomes 0] : 

H ^^pT-{-gM) (20) 

Here we have used Eq. (fTTJ]) to obtain the inequality on 
the RH side. For = 10 13 GeV and m - 0(1) TeV, it 
can be seen that for all couplings in the range 10~ 6 < g < 
1 , we have Ha < Tj ~ m^/g 2 Mp, which is equivalent 
to one-particle decay of the inflaton. This implies that 
the inflaton decay becomes kinematically allowed only in 
the perturbative regime. 

The total decay rate of the inflaton (when we also ac- 
count for decay to the fermionic partner of x) is given 
by r d = .g 2 77i /87r. For 1(T 6 < g < 1 we have T d > H d . 
This implies that the inflaton will immediately decay as 
soon as kinematics allow 13 . The above Eq. (f2TJ!) im- 
plies m 

H d < (irr 3 - ltr 1 ) TeV (21) 



The inflaton decay into two x (and its fermionic partner) is kine- 
matically forbidden for H > H^. However higher order decays 
to light particles via off-shell \ (and its fermionic partner) will 
be kinematically allowed at all times. The lowest order such pro- 
cess is four-body (perturbative) inflaton decay via two off-shell 
X- This decay channel, in addition to phase space suppression, is 
suppressed by a factor of (m^/ypo) 4 - This results in a very small 
decay rate which will not be important for our discussion. 
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for mo ~ 0(1) TeV. This underlines the fact that, re- 
gardless of how large its coupling to x i s : the inflaton 
will not decay until after the flat direction has started 
oscillating, and even then its decay will be strictly per- 
turbative. The expression in Eq. (|21[) for the inflaton 
lifetime is very robust and practically independent from 
how exactly flat direction oscillations terminate (see the 
discussion below). 

As pointed out in Refs. [3, [3] , when the inflaton com- 
pletely decays, the decay products of MSSM fields do 
not thermalize promptly. The flat direction VEV also 
gives masses to gauge bosons and gauginos which slow 
down various scattering processes, i.e. 2 <-> 2 and 2^3, 
required for a complete thermalization. The Universe un- 
dergoes a bout of quasi-thermal phase where the plasma 
obtains (near) kinetic equilibrium, but not chemical equi- 
librium. Eventually when the flat direction starts oscil- 
lating and its VEV decreases the Universe reaches full 
thermal equilibrium. In a model independent case, the 
reheat temperature is governed by the thermalization 

1 /2 r— i 

rate, i.e. Tr ~ (TthMp) , and can be as low as TeV |2j. 
We note that the VEVs required for kinematical blocking 
of preheating, see Table. 1, are sufficiently large to also 
delay complete thermalization of the Universe [2|. 

In some cases the coherent oscillations of the flat di- 
rection can fragment the homogeneous mode into a non- 
topological solitons [45(, known as Q-balls. In which case 
the flat direction VEV vanishes outside the Q-balls, while 
still being large inside them. Then, since the Q-balls oc- 
cupy a very tiny fraction of the space, the inflaton decay 
will become kinematically allowed as soon as Q-balls are 
formed. In principle flat direction oscillations might also 
decay rapidly via preheating due to their gauge interac- 
tions and initial condition ipo 3> mo 14 . An important 
point is that the time scale for fragmentation of flat di- 
rection oscillations, or their non-perturbative decay via 
preheating, lies within the same order as the RH side of 
Eq. (|2~lj) . Therefore, irrespective of the fate of flat direc- 
tion oscillations, Eq. (|2f |) provides a robust lower bound 
on the inflaton lifetime. 



VIII. CONSEQUENCES OF NO-PREHEATING 

The absence of preheating and a delayed perturbative 
decay of the inflaton has interesting cosmological conse- 
quences which we briefly discuss here. 



This does not happen for flat directions which have a non-zero A- 
term (either from higher-order superpotcntial or Kahler potential 
terms). The A-term in this case triggers out-of-phase oscillations 
of the real and imaginary parts of the fiat direction (with compa- 
rable amplitudes). Then the mass of particles which are coupled 
to the fiat direction will not experience a non-adiabatic variation, 
and hence no preheating |46t [47| . 



A. No gravitino problem 

The foremost consequence of NO preheating is that 
there will be no copious production of dangerous relics 
from scatterings in the aftermath of inflaton decay. First 
of all note that there will be no preheat plasma with large 
occupation numbers fk > I, which could lead to disas- 
trous production of relics from scatterings [l|, [28|. On 
the other hand, since the inflaton decays perturbatively, 
the MSSM fields are scarcely populated Q. The chal- 
lenge is to enhance the number density of particles which 
mainly happens via 2 — > 3 scatterings. As pointed out 
in Refs. [l|, 0] the production of dangerous relics such as 
gravitinos is very much suppressed in the initial plasma. 

Large VEVs of the flat directions do not modify the 
non-thermal production of gravitinos during the coherent 
oscillations of the inflaton [I2|, EH, EH • However this will 
not be a threat for BBN. Helicity ±3/2 are not produced 
copiously at the first point, while ±1/2 fermions thus 
produced are mainly inflatinos, which decay along with 
the inflaton (thus long before BBN) [T3 |. 

Another source is direct production of gravitinos via 
the channel inflaton — > in f latino + gravitino (if kine- 
matically open). If the inflaton decay to matter occurs 
very late, this process would have a large branching ratio 
and could lead to overproduction of gravitinos [H], [4^] . 
However this will not be a problem because delayed infla- 
ton decay (due to kinematical blocking) occurs at see 
Eqs. ()20|21|) . in which case the direct gravitino produc- 
tion from inflaton decay is under control. Similar analysis 
holds for moduli fields also since their interaction rates 
are also suppressed similar to the case of gravitinos. 



B. No non-thermal leptogenesis 

Leptogenesis is a scheme for creating a lepton asymme- 
try, which is then partially converted into baryon asym- 
metry via SM sphalerons, for a review see [50(. The 
prospect for thermal leptogenesis is severely hampered 
within SUSY due to late thermalization [2]. This moti- 
vates the case for non-thermal leptogenesis (for example, 
see (H). 

Note that the RH (s)neutrinos can obtain large masses 
through their coupling to H u and L superfields (for de- 
tails see Section TXI Bp . If flat directions including H u 
and/or L develop a large VEV then the resonant ex- 
citation of the heavy (s)neutrinos is kinematically for- 
bidden. If flat directions which do not include H M 
or L develop a VEV, for instance udd, then the RH 
(s)neutrinos obtain large VEV dependent masses from 
the non-renormalizable superpotential through effective 
Yukawa couplings [32J. In either case the resonant pro- 
duction of (s)neutrinos is unlikely for the same reason as 
we discussed above. This seriously severs the prospect 
for non-thermal leptogenesis from on-shell superheavy 
(s)neutrinos. 
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C. No creation of superheavy WIMPS 

In Ref. Q it was advocated that it is possible to ex- 
cite superheavy a weakly interacting massive particle 
(WIMP) from non-perturbative inflaton decay. In most 
of the cases these WIMPs are either charged under the 
SM gauge group, or are a SM gauge singlet which is cou- 
pled to some gauge sector. In either case their production 
through non-perturbative decay of the inflaton will be 
hampered in SUSY. In the former case a MSSM flat di- 
rection developing a VEV induces a large mass to WIMP 
which can kinematically block its production via preheat- 
ing. The latter case is similar to that of right-handed 
(s)neutrinos discussed in the previous subsection. 



D. No Exit from hybrid inflation 

So far we have considered the effects of flat directions 
on preheating which usually arises in models of chaotic 
inflation. Here we briefly discuss the case for models of 
hybrid inflation. 

The simplest SUSY hybrid inflation model has the fol- 
lowing superpotential, see for instance [|| 

WDy* (* 2 - *g) , (22) 

where 3? is the inflaton superfield and * is a superfield 
which contains an auxiliary scalar field ip. Inflation is 
driven by the false vacuum potential y 2x &Q during which 
<f> undergoes slow-roll motion and ip is stuck at the ori- 
gin. Inflation ends when ((f)) reaches the critical value 
4>c = ^o/ (V%y) , at which point the auxiliary field mass- 
squared becomes negative and tachyonic preheating takes 
place [HI . Eventually the two fields settle at (<p) = and 
W = *o- 

Hybrid inflation is often quoted as the most successful 
inflationary model motivated by particle physics, for a 
review see p|. However in all cases the inflaton remains 
a SM gauge singlet. For a successful phenomenology it is 
then imperative that the vacuum energy density during 
inflation is converted into MSSM particles. 

In the simple case given by Eq. (f22|) both the inflaton 
and the auxiliary field are SM singlets. Then the * su- 
perfield can have superpotential couplings to the gauge- 
invariant combinations of the MSSM given in Eqs. (13I4[) . 
Since (ip) = during inflation, any field coupled to ip is 
massless during inflation. This implies that any MSSM 
flat direction <p which includes a field coupled to ip can 
acquire a large VEV, again denoted by ifo, in the infla- 
tionary epoch. Such a large VEV will induce a mass Xipo 
to ip through the following term in the scalar potential: 

A 2 MV- (23) 

If Xipo > y 1 / 2 vl> , the mass-squared of ip will remain pos- 
itive even for (<p) < <p c . As a result (ip) = while <p 
is still slow rolling, and hence tachyonic preheating will 



never take place. More importantly it implies that infla- 
tion will never end! Therefore, unlike chaotic inflation, 
the absence of preheating in hybrid models has a neg- 
ative consequence. Indeed for a graceful exit from the 
inflationary phase one needs to have AtySo < y 1 ^ 2 '$o- 

In more realistic models of hybrid inflation the auxil- 
iary field ip is associated with a Higgs field which spon- 
taneously breaks part(s) of a larger gauge group which 
contains the SM, such as grand unified theory (GUT) 
or an intermediate scale U(1)b-l, upon developing a 
VEV after inflation ends (see for instance [H, [5l|). This 
Higgs field is naturally coupled to some of the fields which 
carry gauge quantum numbers under the larger symme- 
try. Note that this symmetry is unbroken during infla- 
tion since (ip) = 0. Then flat directions which include 
the fields coupled to ip can obtain a large VEV during 
inflation and induce a large mass to it. Again this can 
prevent an end to inflation. 

Therefore graceful exit from hybrid inflation typically 
leads to constraints on the couplings of the auxiliary 
field to gauge non-singlets. Note that such couplings are 
present in any model which is embedded into a parti- 
cle physics model. This must be taken into account in 
any realistic model of hybrid inflation. We will deal with 
some of these issues in a separate publication. 



IX. CONCLUSION 

We argued in this paper that non-perturbative decay 
of the inflaton via resonant or tachyonic instabilities is 
unlikely within SUSY, due to the presence of flat direc- 
tions, and it rather prefers to decay perturbatively. A 
flat direction VEV > 10 13 GeV will be sufficient to kine- 
matically block preheating. 

The key observation is the presence of flat directions. 
Within MSSM there are nearly 300 flat directions, it is 
expected that a number (if not all) of independent di- 
rections would develop large VEVs during inflation. The 
flat directions have Yukawa couplings to the inflaton de- 
cay products and induce large SUSY preserving VEV- 
dependent masses to them. For reasonable (sub Planck- 
ian) VEVs, depicted in Table. 1, this leads to a kinemat- 
ical blocking of the inflaton decay via preheating, even 
for 0(1) inflaton couplings to other fields. 

The inflaton decay will be kinematically allowed only 
after the flat direction starts oscillating, and once its 
VEV has been sufficiently redshiftcd. We found that 
the final decay of the inflaton is perturbative and there 
exists a robust upper bound on the inflaton lifetime 
- (10 - 10 3 )TteV -1 . 

The absence of a violent stage of non-perturbative 
inflaton decay has important implications for particle 
cosmology. Most notably the initial plasma has much 
smaller occupation numbers, i.e. /j < 1, which implies 
that overproduction of dang erous relics through scatter- 
ings is unlikely, see [l|, [2j, |28| . As we have already noticed 
in Ref. Q, the flat direction VEVs suppress the rate of 
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thermalization, therefore the reheat temperature can be 
as low as O(TeV). In which case thermal production of 
gravitinos and other relics will be negligible. 

As a cursory remark we pointed out that SUSY hybrid 
inflation requires a more careful treatment while taking 
into account of the flat directions. Our brief discussion 
suggests that large VEVs of the flat directions can even 
prevent a graceful exit from inflation in these models. 

Finally, since our focus was on the m^cj) 2 case, one 
might worry about preheating in A^> 4 model. First of all, 
in the pure Xcj) 4 case preheating is inefficient and only 
about 5% of the energy density in inflaton oscillations 
is converted into its own quanta [5t| . In any case the 
(j> quanta must eventually decay to MSSM fields, which 
will be kinematically forbidden in the presence of a large 
flat direction VEV as we have discussed. Moreover the 
\<j) A model is rather unattractive as it can produce large 
non-Gaussianity [191 ]. 
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XI. APPENDIX 

A. MSSM flat directions 

There exist a large number of directions in the field 
space of supersymmetric theories, known as flat direc- 
tions, along which the scalar potential identically van- 
ishes in the limit of exact SUSY. In this limit the scalar 
potential of MSSM, denoted by Vmssm, is the sum of the 
F- and D-terms and reads 

V = Y.\ F i\ 2 + lY,9 2 a DaDa > (24) 

i a 

where 

^^SM D » =x * T a x . . (25) 

Here the scalar fields, denoted by x^, transform under a 
gauge group G with the generators of the Lie algebra and 
gauge coupling are by T a and g a respectively. 



For a general supersymmetric model with TV chiral su- 
perfields, it is possible to find out the directions along 
which the potential in Eq. (fM)) vanishes identically by 
solving simultaneously 



Xi I'm Xj , 



dW 



0. 



(26) 



Field configurations obeying Eq. (|26|) are called respec- 
tively D-flat and F-flat. 

D-flat directions are parameterized by gauge-invariant 
monomials of the chiral superfields. A powerful tool for 
fin ding the flat directions has been developed in [H, [H, 
[H, I591I60I. IfJH , where the correspondence between gauge 
invariance and flat directions has been employed. In 
particular all flat directions have been classified within 
MSSM [H|. 

Adding the inflaton superfield which is a SM singlet 
does not affect D-flatness. However one might worry that 
the inflaton coupling to matter would ruin the F-flatness 
as new terms can now appear in the superpotential. For 
example consider the non-renormalizable superpotential 
terms: 



— *H u Qu, — *H d Qd, — *H d Le 

M P ' M P ' M P 

— *QLd, — *udd, — *LLe. 

M P Af P M P 



(27) 



which can arise in addition to the renormalizable one in 
Eq. ([5]). If W = H d , see Eq. ©, terms in the first row 
preserve i?-parity while those in the second row violate 
it. The reverse situation happens if = L. If i?-parity is 
a discrete subgroup of a gauge symmetry, it will remain 
unbroken by gravitational effects. In this case only those 
terms in Eq. (|27[) which preserve i?-parity can appear in 
the superpotential. On the other hand, like other global 
symmetries, i?-parity will be supposedly broken due to 
gravitational effects if it is not protected by some gauge 
symmetry (see for example Ref. [57|). In this case the 
superpotential can include all terms in Eq. (|27|) regardless 
of what ^ represents. 

Further note that in this case the LSP decays via terms 
which violate i?-parity. However the decay is suppressed 
by both Mp and m^, and hence the LSP lifetime is much 
longer than the age of the Universe. These interactions 
can also led to proton decay though at time scale which 
is much longer than the experimental bound. Therefore 
i?-parity breaking through Mp suppressed operators will 
not be constrained. 

In models of large field inflation, in which (4>) > Mp 
during inflation, the terms in Eq. (127)) result in effective 
renormalizable superpotential terms. Note that terms 
from the first row are exactly the same as the SM Yukawa 
couplings, and hence do not lead to any new constraints. 
On the other hand, the second row results in terms which 
are absent in the MSSM superpotential. Therefore, if 
allowed, they will lead to new F-flatness constraints and 
lift some of the flat directions. 
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Note however that a large subset of MSSM flat direc- 
tions survive. In particular those which are only made 
up of Q, u, e superfields will not be affected at all . Inter- 
estingly these directions are coupled through the MSSM 
superpotential in Eq. @ to both of the H u and \& super- 
fields, whether *f? = or = L, see Eq. ([6]). As we will 
see, this is very important for shutting off preheating. 

B. Inflaton couplings to SM gauge singlets 

Inflaton being a gauge singlet need not directly couple 
to the (MS)SM sector, but instead can do so through an- 
other SM gauge singlet. Phenomenologically the best 
motivated example is when the inflaton is coupled to 
the RH (s)neutrinos. Note that the couplings of RH 
(s)neutrinos to the (MS)SM sector can explain the origin 
of light neutrino masses via see-saw mechanism |52j | . 

Assuming that the RH (s)neutrinos obtain masses from 
some other source, the relevant part of the superpotential 
will be given by: 

W D ^m ** + ig*NN + /iNH„L + ^M W NN . (28) 

Here <f>, N, L, H„ stand for the inflaton, the RH neutrino, 
the lepton doublet, and the Higgs (which gives mass to 
the top quark) superfields, respectively. Also, and 
M/v denote inflaton and RH (s)neutrino masses, respec- 
tively. Here h denotes the 3x3 neutrino Yukawa matrix. 
For simplicity, we have omitted all indices in h matrix 
and superfields, and work in the basis where Mjy is diag- 
onal. The inflaton coupling to the RH (s)neutrinos can 
be quite large: g ~ 0(1). 

Any flat direction that includes H u and/or L can in- 
duce a large mass to N through renormalizable couplings. 
Moreover note that N can also couple to MSSM fields 
via non-renormalizable superpotential terms the same as 
Eq. (|27p with <fr replaced by N. In the presence of large 
flat direction VEVs, non-renormalizable interactions in 
Eq. (|2"T|) can lead to large effective couplings between 
N and MSSM fields \3^. For a reasonable flat direc- 
tion VEV, one can have A c ff > 3 x 10~ 4 . Then non- 
perturbative decay of the inflaton to RH (s)neutrinos will 
be kinematically forbidden similar to our earlier analysis 
as in Section [TVl 

C. Categorizing flat directions 

A close inspection to the MSSM flat directions shows 
that they belong to one of the following groups (for in- 
stance see [24|]): 

1- Directions which include two or more Q from dif- 
ferent generations: QQQL. These directions have 
Yukawa couplings > 10~ 3 to both H„ and both. 

2- Directions which includes two or more u from dif- 
ferent generations: uude and uuuee. These direc- 
tions have Yukawa couplings > 10 -3 to H u . 



3- Directions which include two or more d from dif- 
ferent generations: udd and LLddd . These direc- 
tions have Yukawa couplings > 10 -3 to H^. 

4- Directions which include Q and u from different 
generations: QuLe and QuQue . These directions 
have a coupling > 1CP 3 to H„. 

5- Directions which include Q and d from different 
generations: QdL. These directions have Yukawa 
couplings > 1CP 3 to H<j. 

6- Directions which include two L from different gen- 
erations: LLe. These directions have Yukawa cou- 
plings > 1(T 3 to H d . 

Note that the requirement that more than one genera- 
tion of squarks and/or sleptons be involved comes as a 
direct consequence of F- and D-flatness [22j ■ Now let us 
consider the SI/ = and VP = L cases separately, see 
Eqs. 460). 

* = ELj. In this case Ai > 1CT 3 and/or A 2 > 10~ 3 for 
all of the MSSM flat directions, see Eq. (fTO]) . implying 
that A > 3 x 1CT 4 . 

* = L. In this case Ai > 1CT 3 , hence A > 3 x 1(T 4 , 
for flat directions listed in 1, 2, 4. There are exceptional 
flat directions for which A <C 1CP 3 : 

• udd: there are three such directions (with u be- 
longing to the first generation) for which Ai ~ 
0(1(T 5 ) and A 2 = 0, hence A - 0(1CT 5 ). 

• QdL: there are 6 such directions (with Q belonging 
to the first generation) for which Ai ~ 0(1O -5 and 
A 2 = 0, hence A - 0(1CT 5 ). 

• LLddd: there is one such direction (where the two 
L are orthogonal to St. As mentioned earlier, tp 
cannot acquire a large VEV since its coupling to 
the inflaton induces a mass 3> Hi.) for which Ai = 
A 2 = 0, hence A = 0. 

• LLe: there is one such direction (where the two L 
are orthogonal to \P) for which Ai = and A 2 ~ 
0(1(T 5 ), hence A - 0(1(T 5 ). 

The exceptional directions constitute a small subset of 
all MSSM flat directions: 11 out of nearly 300. More- 
over, considering multi-dimensionality of the space of all 
flat directions, it is very unlikely that a VEV grows ex- 
actly along one of the exceptional directions during in- 
flation. To elucidate consider flat directions represented 
by the udd monomial. When different generations are 
taken into account, this monomial represents a space of 
complex dimension 9 [12]. While, exceptional directions 
span a three-dimensional subspace. Therefore even with 
a probabilistic argument it is very hard to imagine that 
a non-zero VEV will be confined to this subspace of flat 
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directions. For the bulk of 9-dimensional space the super- implying that Ai > 10 3 and A > 3 x 10 
field, u, has comparable components from all generations, 
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